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The critical lines of binary mixtures of conformal fluids are calculated on the basis of the one-fluid Van 
der Waals theory. The conditions for the existence of the gas/gas equilibrium of the first and second type are 
shown and discussed. The analytical expressions for partial derivatives of Gibbs free energy with respect to 
concentration, up to the third order, are given in the appendix. 
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1. Introduction 

In this work we study the binary mixture at the 
gas/liquid critical line (the Plait point curve). Our emphasis 
will be on locating the line rather than discussing the 
behavior of the critical point itself [1] : it is felt that the 
former problem is of scientific and practical interest which 
has not yet been fully explored. Our procedure is based on 
one fluid Van der Waals conformal solution model since it 
is now well established [2-6] that phase equilibria in 
general can be predicted successfully by this method. 

The organization of the paper is as follows: the 
corresponding states approach is recalled briefly in section 
2 and the calculation of the critical lines is discussed in 
section 3 with the lengthy working equations relegated to 
the appendix. The phenomena of gas/gas equilibria are 
next discussed. In particular, some variations of the critical 
line over a range of pressure, P 9 temperature, 7, and mole 
fraction, x, for various types of mixtures are described in 
section 4. Brief comparisons with the more restricted 
approach from simple equations of state [7,8] are given. 

Throughout this discussion, argon is used as the 
reference substance [9]; results with other reference 
equations will be published in a forthcoming paper [10]. 

2. Conformal Solution Theory 

A basic postulate of the Van der Waals one fluid theory 
is that if the components of a mixture — characterized by 



the usual parameters, € and cr — separately obey 
corresponding states, then their mixture will also obey 
corresponding states as if it were a single substance. This 
hypothetical equivalent substance is characterized by some 
suitable composition dependent average parameters, (e) 
and (cr) . The statistical mechanical basis behind the 
theory has been discussed by several authors, in particular 
by Henderson and Leonard [11], Leland, et al. [12], Jacucci 
and McDonald [13], and by Hanley and Evans [14]. The 
concept cannot be strictly valid but the model is 
nevertheless appealing and often successful in practice. 

Accepting that the singie component equivalent 
substance is a valid entity, its properties are then 
determined from those of a known reference substance via 
corresponding states. See, for example, references 3 and 5. 
For the Gibbs free energy the expressions are 



GJTJ>jc) = G X (T,P) + RTlx^nx, 

i 

G X (T,P) = f£ {T/f x , Ph x /f x ) - RTAh x 



(1) 

(2) 



Here the subscript m refers to the mixture, x to the 
equivalent substance and o to the reference substance; 



f. = 






and A = 



<<r-' ! > 



'Su,, |)( ,rted by .he Resean-h Projec. 03.10.1.02.01 of .he r<,ii s h Academy of are p arame ter 8 of the equivalent substance with respect to 
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h x = 2 xiXjhij 
f x = (lxiXjfijhij)/h x 



(3a) 



(3b) 



where 






f it 



fa = 






1 n 



The parameters referring to unlike pairs are usually 
expressed in terms of a mixing rule: 



h v = ^v 



|*a+i*, 



(3c) 



ft = iAfufjJ 

where 7]- and £•■ are semi-empirical correction factors which 
are close to unity; V& is the critical volume and T[ t is the 
critical temperature of the component i. 

3. Critical Lines 

Phase behavior of a mixture can be characterized by the 
course of the gas/liquid critical line in (P—T—x) space 
joining (or not) the critical points of the pure components 
[15]. The line depends on the interactions between unlike 
(i-j) molecules and on the relationship between the critical 
temperatures and critical volumes of the pure components. 
Many investigations have examined this question both from 
the theoretical and experimental viewpoint. See, for 
example, references [16] and [17]. Of special interest, Van 
der Waals [18] and Kamerlingh Onnes and Keesom [19] 
predicted that a phase separation could occur in a binary 
mixture at temperatures well above the critical temperature 
of either component. Hence the term gas/gas equilibria. 
[The effect was verified experimentally in 1940 for the 
system NH 3 -H 2 [20] and one of the earlier theoretical 
discussions is reported in reference [21] for the system Xe- 
He.] 

The wide range of possible phase behavior is covered by 
the extensive work of Scott and Van Konynenburg 
[7,22,23] who used the Van der Waals equation of state 
with the mixing rules 

<*m = a U x2 \ + 022*2 + 2^*2^2 
b m = ^11*1 H~ ^22^2 + ^ X ] X 2 b u 



where a- and b- are the Van der Waals constants, to 
calculate critical lines as functions of 7" c , V c or of a and b. 
Further calculations of critical lines based on the Van der 
Waals and Guggenheim reduced equations of state were 
later reported and compared with experiment [8,24,25]. 

3.1. Equations 

The critical point of a binary mixture at a given T,P and 
x is defined by the equations of diffusional or material 
stability [26]: 



U 2 / 

\ a * 4 A 



= o 



= 



> o 



(4a) 



(4b) 



(4c) 



where Y m = GJRT, with G m , the Gibbs free energy of the 
mixture: x(=x 2 ) is the mole fraction of the second 
component. The locus of the critical points is of course the 
critical line. Note that there is no formal difference between 
a liquid/liquid curve and a gas/liquid (plait point) curve. 

Equation (4) can be solved in two ways, the first 
graphically. The Gibbs free energy is plotted versus mole 
fraction at T and P: the common tangent determines the 
boundary of the phase separation and a critical point is the 
T,Pjc point at which the common tangent vanishes. The 
second way is to solve equations (4a) and (4b) directly, with 
the constraint of equation (4c). This is the procedure used 
here, coupled with the one fluid model. 

Accordingly equations (1) and (2) are expressed in the 
following form: 

YJTfj;) = Y (T o ,P) - fnh x +xtnx+ {l-x)fn(l - x) 



Y m (TJ>*) = GJTJ>*)IRT (5) 

Y„(T A) = G {T o ,P o )/RT o , T o = T/& P = Phjf x 

Equations (4a) and (4b), which are to be solved, are 
obtained explicitly by differentiation of equation (5) with 
respect to x — the mole fraction of the second component 
resulting in, 
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Making use of equations (4a) and (4b) we obtain a 

v t'2 i 9v „p' t' . v d'2 , v t" necessary condition to satisfy equation (8). Namely at the 
Yo,T,t1 T £ioj> y Tr l o -f Y j>j>r -f Yo,Tl . . , . y ^ ^ 

critical point 



T,P 



+ r ,/»;/ _ a;//z x + (h'Jhf + i/*(i -*) 



(6a) 




= r ,T.T,rT' * + 3K,.r,^/ > ;r; 2 + 3r oW />; 2 r; 



r^ 



+ 1W,pPJ 3 + SYomT'oT: + 3Yo,tAPoT2 + /W) 



<^ mr =0 



(9a) 



Km.wdr + Ymj^rdT + Ymm>dP = (9b) 

/77 r . 



(9c) 



+ SYoj>j>PoPo + Y .tTS + WV - Ai'VA, (6b) E( I uation < 9c ) follows from equation (9b) because we are on 

the critical line, equation (4b). There is no reason for Y mxaJ 
to be infinite, thus to satisfy equation (8) we must have 

+ Sh' s h"fh 2 x -2(h'Jh x f- l/x 2 + 1/(1 -xf 



Y„„ P = 
where Y oTTi etc., are defined in the appendix, the prime 

denotes the derivative with respect to x. The equations were 

solved numerically using the Newton method. u f«ii^*„« *u„* ;r *u * ™ c 

J & It tollows that it the system of equations: 

3.2 Gas/Gas Equilibria 



(10) 



Equation (4) and the one-fluid modifications are quite 
general and we have applied them to several systems of 
interest; for example, to mixtures of the noble gases, to the 
mixture N 2 /CH 4 and to the mixture N 2 /C 2 H 6 . These results 
will be discussed in a later publication [2]. Here, however, 
we report on the phenomenon of gas/gas equilibria for 
model systems. We refer to the classification of Scott [7]. 
Type I: Equilibria are classified by the behavior of the 
critical line in the P-T projection. For type I behavior, the 
line must begin from the critical point of the less volatile 
component with a positive slope, i.e., 



(7) 



where subscript c denotes differentiation along the critical 

line. 

Type IT. For this class, the P-T projection has a minimum 

with respect to temperature as the pressure increases. 

Therefore a condition is 



Y = 

m.x+x,x 



y ., = o 



(id 




(8) 



has a solution at a temperature below the critical 
temperature of the less volatile component, then gas/gas 
equilibrium of the second type occurs. 

In order to calculate Y mfXfX j>, we derive from the one-fluid 
theory: 



= Y ,T,Tj>T h x /f x + 2Y Of P,TfP' Tfa'f* 

+ YosAZP'oK/P + T» h*/f x ) + Yoj>j>j>Po 2 hJfi 

(12) 

+ y ^p*ip + P'X/f x ) + y <kP p:;/p 



Mixed Type: A third possibility can easily occur from our 
calculations, even though it has not been observed 
experimentally; namely a mix of type I and type II is 
permissible. For instance the derivative (dP/dT) ( . can be 
positive initially but then a, maximum can occur with 
respect to the temperature at a relatively low pressure, and 
afterwards a temperature minimum is found at higher 
pressures. 
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4. Results 

We calculated a critical point at a given state point via 
equations (4)-(6) — using the mixing rules of equation (3) 
with r\ lJ and i; lJ set equal to unity — so that equations (7) 
and /or (8) were satisfied. In fact equation (6) was solved for 
a given fixed ratio of the critical volume of the pure 
components, and the ratio of the critical temperatures then 
varied until the minimum value of this temperature ratio 
was found which would satisfy equation (7) or (8). Table 1 
lists the results, i.e., the appropriate value of T 22 /T c u for 

fixed values of V C 22 IV C U between 0.7 and 4.5. As a matter of 
interest the overlapping results of Hicks and Young [24] are 
included. 

TABLE I. Minimal value of T^/In for gas/gas equilibrium. 



1 22/1 11 


First kind 


Mixed kind 


Second kind 




(Type I) 




(Type II) 


0.7 






1.88 


0.8 






2.07 


0.9 






2.25 


1.0 


7.06 


4.8 


2.42 




[14.0 ] a 




[2.5 r 




[81.0 f 




[4.0 ] b 


1.05 


7.40 






1.1 


7.79 






1.15 


8.22 






1.2 


8.73 






1.25 


9.38 






1.3 


10.4 






1.35 


14.2 






1.5 




5.43 




1.66 






3.47 


2.0 




5.86 


3.95 


2.5 




6.18 


4.60 


3.75 






5.85 


4.0 




6.73 


6.06 


4.2 






6.23 


4.5 




6.87 





a Calculated [24] by means of the Guggenheim equation of state. 
b Calculated [24] by means of the Van der Waals equation of state. 



The results are perhaps clearer when shown in figure 1. 
Curve I is the plot of the temperature ratio in the second 
column of table 1 versus V 22 IV c n . It represents a limiting 
critical point line above which and including which 
unambiguous first kind, or type I, of gas/gas equilibria 
occur. In other words if the temperature and volume ratios 
of the components of a binary mixture are such that they 
lie inside curve I, the Van der Waals one-fluid model will 
predict gas/gas equilibria of type I. Curve II is the limiting 
line above which and including which equilibria of type II 



T C if 

i 22 / \ }] 




V c /v c 

v 22 ' V 11 



FIGURE 1. Ratios of the temperature at fixed ratios of the critical volumes 
designating regions of different classes of behavior of the critical lines for 
gas/gas equilibria. 

Curve I — the lower limit of appearance of the gas/gas equilibrium of the first kind or type I [7]. 
Curve I A — the lower limit of appearance of the gas /gas equilibrium of the mixed kind. 
Curve II — the lower limit of appearance of the gas/gas equilibrium of the second kind 
or type II (7J. 



will occur. Curve IA represents equilibria of the mixed 
kind. This curve is somewhat arbitrary since we required 
that equation (7) be satisfied at pressures between P 22 and 
P 22 + AP is arbitrary set at 7.5 MPa, and at temperatures 
greater than the critical of the less volatile component. We 
have determined, however, that curve IA with AP=7.5 
MPa represents the upper boundary: Use of AP<7.5 MPa 
simply tended to collapse IA towards II. 

Figures 2 and 3 show sample plots of the critical lines in 
the P-T projection. Curve I of figure 2 is typical of type I 
behavior: in this case with V c n IV c 22 =\3 and T € n /T c 22 = 
10.4. Curve IA is typical of a mixed kind result with the 
volume and temperature ratio of 2.5 and 6.18, respectively. 
Figure 3 gives type II results for the ratios (1 and 2.42) 
and 1.25 and 2.59. 

5. Discussion and Conclusions 

Our first remark is that we have proposed a very 
efficient general procedure to calculate the critical lines of a 
mixture based on the one-fluid Van der Waals conformal 
solution theory. Computational time is rapid. As an 
example of the procedure we report results for the critical 
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Figure 2. P-T projection of type 1. at Vl 2 /Vli = 1.3 and T§ 2 /Tf 1 = 10.4, 
and of the mixed type. (la) at V§ 2 /Vf2=2.5 and T§g/Tfi = 6. 18. 
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Figure 3. 
TWTh 
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P-T projection for type II behavior at Vfa/Vfi^l.O and 
= 2.42 (curve a), and V$ 2 /V c n = 1.25 and T^ 2 /T c n =2.S9 (curve 
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lines in gas/gas equilibria, but we have used the method to 
estimate gas/liquid and liquid/liquid lines for several 
hydrocarbon systems with success [2,25]. 

We next remark on the relationship of our curves with 
experiment. Unfortunately this is not clearcut for the 
following reasons: 

(a) Conformal solution theory is not strictly valid and 
tends to break down for a mixture whose species differ 
widely in size, mass, and interaction energy [13,14]. 
Unfortunately such mixtures are the more likely to show 
gas/gas equilibria, especially mixtures with helium. Our 
predictions of type I behavior for mixtures of He-C 2 H 6 , He- 
C0 2 and He-Xe (for which VyV c u ~4 and 7^/7^-30) 
where not too satisfactory. Predictions of type II behavior 
for the mixtures Ne-Kr, Ne-Xe and Ne-CH 4 were somewhat 
better. Also the model correctly gave the negative result 
that mixtures whose components are not disparate, e.g., Ar- 
Kr, Ar-Xe and Kr-Xe, do not show gas/gas equilibria. 

(b) Real mixtures do not obey corresponding states, 
although the method can be generalized to include the so- 
called shape factors introduced by Leland and coworkers 
[6]. 



(c) The procedure — common with all mixture procedures 
available at this time — requires interaction parameters such 
as the £ 12 and 7] l2 of equation (3) and the numerical 
properties of a mixture are often, but not necessarily, 
extremely sensitive to a particular choice of £ and 7). We 
have set £ and 7] equal to unity but their influence on the 
critical lines is quite pronounced, which is demonstrated in 
figures 4 and 5. The general tendency is: lowering £ >n shifts 
the critical line towards a higher temperature, (the same 
qualitative effect is obtained by increasing the critical 
temperatures ratio with the critical volumes ratio held 
constant): increasing T} 12 increase the curvature of the 
critical line. 
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Bowie for preparing the manuscript. 
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Figure 4. Effect of^ 12 and 7] 12 on the projected P-T curve of the mixed kind. 
All calculations at ^22/^11 = 2. 5 and 7 22^ 11 = 6.18 (see fig. 2). Solid 
curve: ^ 12 = T) 12 = 1.0, dashed curve: % 12 =1.0, t] 12 = 108, dotted curve: 
Vi2=1.0, kn = 0,8. 



Figure 5. Effect of i; 12 and t] 12 on the projected P-T curve of type II. 

All calculations at f / 22 // ^ii=2. 18 and T22^Ti l = 4.71. Solid curve: 
r li2 = Zi2 = 1-0, dashed curve: %i2=l, 7)i2 = ll, dotted curve: r) 12 = 1, 
£ 12 = ().9, dash-dotted curve: 7) 12 =1.1, ^ l2 = 0.9. 
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Y oJ >j = -V/RT 2 + 




IRT, 





Y oJJ = 2H(T,P)/RT 3 



(dS\ h 2 S } 

— \/RT 2 J- 
*tI p \ar/„ 

YoATs == 2V/RT*-2 ( — /RT 2 +[ | / Rl 

dT L \ar 2 



*+w = 



*oJ> % PJ> 




/RT 2 



r d 2 V 



dPdT 



/RT 



where we have omitted subscript o on the right hand side 

for brevity. 

The equation of state of argon is: 



P = P(T,p) 



l/V 



(Al) 



7. Appendix 

The partial derivatives of the function Y (T o ,P ) which 
are necessary for calculations of the second and third 
derivative of Y m (T,Pjc) with respect to x are given in this 
appendix. H is the enthalpy, F the Helmholtz free energy 
and S the entropy. We define Y oljk as the "derivative of Y 
with respect to /, j and k, 

Y oT = - H(T,P)/RT, 

Y oJ> = V/RT, 



In this work P and T are used as the independent 
variables, thus for any given pair, equation (Al) is solved 
for p o . 

Calculation of H ( T ,P ( ). 

GJ e XT,P) = F r **(T,V) + RT(Z- l)-RTSnZ 

Z o = PVIRT (A2) 

From = j p o[P(T,pvp' 2 -RT/ P ']d P ' 
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After integration of equation (Al) accordingly to equation 
(A2) we obtain Fo es (T,V) as a polynomial in To and p . The 
superscript "res" denotes the residual functions defined in 
general by X t e *=X(TJ 3 )-X id (T 9 P): 



$ "(T,V) = 



J dF res (T,V)^ 



dT 



y; s (r,P) = 5 res (r,F) + RTfnZ 



G fl (T,P) = G Tt \T,P) + RT/nP/kT 



S (TJ>) = STiTJ 3 ) - R( fnP/kT - 1) 



(A3) 



(d 2 S (T ,Vj\ 



(* 2 sr{T ,v o \ 



this quantity is obtained by differentiation of equation (A3). 
Calculation of other derivatives. All derivatives of the 
pressure, P o9 are obtained by direct differentiation of the 
equation of state (Al). The derivatives of the volume V 
with respect to the temperature and pressure are 
reexpressed by the derivatives of the pressure P : 



H£TJ>) = G res (T,P) + TS res (T,P) + RT 



(A4) 



Calculations of 



fd 2 S (T o ,P^ 
and 



dT„ 



a 7? 




* 2 V. 



LaP? 



teJLTW 



fdsrr(T ,ry 






(A7) 



(A8) 



(A9) 




dT. 



e7\ 



this quantity is obtained by direct differentiation of the 
equation (A3) with respect to the temperature. 



+ 2 



r l \dTj \dT 2 1 






(A10) 



(All) 
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